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Abstract. The secular instability and nonlinear evolution of the m = 2 f-mode (bar- 
mode) driven by gravitational radiation reaction in a rapidly rotating, newly formed 
neutron star are reviewed. There are two types of rotating bars which generate quite 
different gravitational waveforms: those with large internal rotation relative to the 
bar figure (Dedekind-like bars) have GW frequency sweeping downward during the 
■ evolution, and those with small internal rotation (Jacobi-like) have GW frequency 

sweeping upward. Various sources of viscosity (which affects the instability) in hot 
nuclear matter are reexamined, and the possible effect of star-disk coupling on the 
bar-mode instability is also discussed. 

o 

I INTRODUCTION 

6, 

Nonaxisymmetric instabilities can develop in rapidly rotating fluid bodies when 
the ratio (3 = T/|VK| of the rotational energy T to the gravitational potential energy 
W is sufficiently large (e.g., Chandrasekhar 1969). In particular, the I = m = 2 
f-mode (Kelvin mode), or bar- mode, becomes dynamically unstable when f3 > 
(3d yn — 0.27. This /3d yn , originally derived for incompressible Maclaurin spheroid, 
is relatively insensitive to the equation of state and differential rotation (Pickett et 
al. 1996; Toman et al. 1998), although it tends to be reduced by general relativity 
(Shibata et al. 2000). The consequence of the dynamical bar-mode instability has 
been extensively studied using numerical simulations: the mode grows to nonlinear 
amplitude by shedding mass and angular momentum from the ends of the bar in 
the form of two-armed spiral pattern, and the central star assumes a bar shape that 
lasts many rotation periods (e.g., Tohline et al. 1985; New et al. 2000; Brown 2000; 
Shibata et al. 2000). It is important to note, however, that these simulations start 
out with a stationary, dynamically unstable star and such an initial condition may 
not be realized in an actual core collapse (e.g., Rampp, Miiller and Ruffert 1998). 

However, it is very likely that in a rotating core collapse (or neutron star binary 
merger), after the messy dynamics is completed, the newly formed NS settles down 
into a dynamically stable configuration, with (3 < 0.27, and yet still suffers further 



secular, rotational instability. "Secularly unstable" means the instability is slow, 
and is driven by gravitational radiation reaction. In this paper, we will focus this 
secular instability and the possible gravitational wave signals from newly formed 
neutron stars in the first seconds/minutes of their lives. 



II CFS INSTABILITY: F-, G- AND R-MODES 

The gravitational radiation driven instability, or Chandrasekhar-Friedman- 
Schutz (CFS) instability (see e.g., Friedman 1998 for a review), arises for the 
following simple reason: Suppose in the rotating frame of the star we set up a 
perturbation (mode) (with frequency u r as in e im <t>+ lu) rt^ travels opposite to 

the rotation. In the inertial frame, the perturbation will be dragged backward by 
the rotation, and the mode frequency becomes u>i = uo r — mQ s . If the spin Q s is 
sufficiently large, the perturbation will be prograde in the inertial frame. Since the 
mode has negative angular momentum (because the perturbed fluid does not rotate 
as fast as it did without the perturbation), as the mode radiates positive J through 
gravitational radiation, the mode's angular momentum will be more negative; that 
means the mode is unstable. 

The CFS instability mechanism can be applied to different modes of NSs: 
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FIGURE 1. Growth times of gravitational radiation driven CFS instability for different modes 
in a neutron star as a function of (3. 



1. F-modes: The f- modes are fundamental acoustic waves, corresponding to a 
global distortion of the star. For the m = 2 mode (bar-mode), the instability occurs 
when (3 > /3 SCC = 0.14 — This critical (3 SCC is only slightly affected by equation of 
state, and is somewhat reduced (to as small as 0.1) by strong differential rotation 
(Imamura et al. 1995) and by general relativity (Stergioulas & Friedman 1998). 

2. R-modes: In the last few years it was realized that the so-called r-modes are 
always secularly unstable (for stars consisting of inviscid fluid) (Andersson 1998) 
and the growth time of the mode through current quadrupole gravitational radia- 
tion can be interestingly short (Lindblom et al. 1998). A lot of works are currently 
being done on r-modes (see http:/ /online. itp. ucsb.edu/online/neustarsOO/si-rmode- 
sched.html for a recent meeting devoted to r-modes), but they are beyond the scope 
of this paper (see Ushomirsky's paper in this proceedings). 

3. G-modes: G-modes are also unstable when Q s is (approximately) greater 
than the mode frequency (of nonrotating stars), which typically is ~ 100 Hz, or 
0.1f2 max . However, the growth time (despite the finite quadrupole moment) is quite 
long, an thus not interesting (Lai 1999). 

Figure 1 shows a comparison of the CFS instability growth times for different 
modes in a M = 1.4M , R = 20 km neutron star. We see that for sufficiently 
large (3, the m = 2 f-mode (the bar-mode) is the most unstable mode, with the 
shortest growth time (of order seconds to minutes). So the bar-mode instability 
is most robust if a neutron star is formed with /3 in the range of 0.2 to 0.27 - 
numerical simulations indicate that this is indeed possible (Zwerger & Miiller 1997; 
Rampp et al. 1998; see contributions by New and Brown). In the remainder of the 
paper we will be concerned with the bar-mode only. 

The growth time shown in Fig. 1 is due to gravitational radiation only. There 
are several complications that can affect the net growth time. We will discuss two 
of them in the next two sections. Other issues, including the effect of magnetic 
field, can be found in Ho & Lai (2000). 



The first complication concerns the viscous dissipation, which tends to suppress 
the GR driven instability. For a proto-NS with T <; 1 MeV, the shear (kinematic) 
viscosity due to neutron-neutron scattering, v = rj/p ~ 14p^ 4 T^g V cm 2 s _1 (Flow- 
ers & Itoh 1976), is negligible (since t visc ~ R 2 /u is much greater than the growth 
time of the mode, £gr)- The shear viscosity due to neutrino- nucleon scattering is 
i] = n v pJ u /5 (Goodwin & Pethick 1982), where n u ,p u ,l u are the neutrino number 
density, momentum and scattering mean free path. Thus 



where we have used l v ~ 2000(^/30 MeV) 3 cm. Contrary to some earlier claims 
(e.g., Lindblom & Detweiler 1979), the neutrino shear viscosity is also negligible. 



Ill VISCOSITY IN PROTO-NEUTRON STARS 




(1) 



The most relevant viscosity in a proto-NS is the bulk viscosity. This arises 
from the fact that in an oscillation which involves compression, the matter will be 
temporarily out of chemical equilibrium (in this case, /^-equilibrium between n,p,e) 
and it will try to relax back to the equilibrium by the weak processes 

e + p^n + v e , n -> p + e + u e , (URCA) (2) 

or 

e + p + N->n + N + v e , n + N-^p + N + e + u e , (Modified URCA) (3) 

and therefore emitting neutrinos; this neutrino emission serves as damping of the 
oscillation. 

The standard bulk viscosity widely used in the last decade has been the one 
derived by Sawyer (1989) (correcting a factor of 100 typographic error) 

C = 1.5 x 10 32 p 2 15 T^ cY uj- 2 g/(cms), (Modified URCA) (4) 

where T Me v is the temperature in MeV, p = 10 15 pi 5 g/cm 3 is the density, and u is 
the mode (angular) frequency in s _1 . This is quite large, with the corresponding 
damping time R 2 p/( ~ lOT^fy s for oj ~ 10 3 s^ 1 . Thus one may conclude that 
the bar-mode instability is suppressed for Tyiev *S 1 ( see Ipser & Lindblom 1991 for 
more detailed computations). Similarly, if direct URCA process operates, equation 
(4) should be replaced by 

C ~ 1.5 x 10 36 p^McV^" 2 g/(cms) (URCA) (5) 

where we have used the free npe gas relation n e /no = 0.0765 p\ b (no = 0.16 fm~ 3 ) to 
evaluate the equilibrium electron density (this is not consistent since in a free npe 
gas URCA precesses are suppressed; but it is adequate for estimate). Equation (5) 
is valid for cu ^> 118x _2//3 T^ eV s -1 (where the electron fraction a; = n e /n = 0.02p 15 ). 

However, one should be careful when using eqs. (4)- (5) at high temperatures. 
In fact, when T becomes sufficiently large, the bulk viscosity must go down. The 
reason is that as T increases, the timescale to relax back to /9-equilibrium be- 
comes shorter than the oscillation period. So the matter will stay very close to 
/9-equilibrium during the oscillation and there is very little extra neutrino emission 
associated with the oscillation. Another effect that needs to be included at high 
temperatures is neutrino absorption (such as v e + n — > e + p) which also helps to 
speed up relaxation to /9-equilibrium. 

We now outline the derivation of the neutrino bulk viscosity in the regime where 
matter is opaque to neutrinos (see Sawyer 1980; Lai 2001). Consider the emission 
and absorption of z/ e 's. The u e distribution function / satisfies the Boltzmann 
equation 



(6) 



where j is the emissivity, and K — jll + e T J is the absorption cross section per 
unit volume (corrected for the effect of stimulated absorption), 5/i = \i e + /i p — /i n . 
In equilibrium (5/j, = 0; we assume [i v — 0), we have f — fo — feq — (l + e ^ ) 
and jo = K f . In a perturbation with 5f oc e tuJt , we find 



(iu; + K )8f = 



(1 + e E / T ) 2 \ T 



5fi\ 
7^ K o- 



(7) 



Similar consideration for v e gives (iuj+K,o)Sf = —e E l T {\ + e E ^ T ) 2 (S/i/T) k , where 
we have used jo = jo and Rq = kq. The variation of electron fraction x = n e /n is 



A e ff, 



(8) 



where A c ff = T 2 k /6. Since the matter is degenerate to a good approximation, 5/i 
depends only on p, x. Thus 
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FIGURE 2. Bulk viscosity ( as a function of temperature. The lower temperature segments 
correspond to (neutrino) optically thin regime, and the high temperature segments to the opaque 
regime. The density is pis = 1, and the mode frequency u/(2ir) = 10 3 Hz. 



where h c q = k + (\ c s/n) (d5fj,/5x) n . The energy dissipation rate per unit volume 
(averaged over oscillation period) can be calculated by 



<*>-J(''£)->>^[<-< 

Using eq. (9) we then find the bulk viscosity 



M?) ' (10) 



A cff / dSfx\ ( dP\ \ cS ( fx n 
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C ^ + Kl s \dn)\dx) p c^ + ^UJ' (U) 

where the second equality applies for free npe gas (fi n = 140 p% 3 MeV is the neutron 
Fermi energy), for which (d5/i/dx) n = 3(n/n e )fi n and thus k b q ~ k . 
For modified URCA process, we have 

«o = 1.034 xlO^T^vS" 1 , (12) 
C = 7.9 x 10 31 p 2 l5 T^ cV {u 2 + g/(cms). (Modified URCA) (13) 

If direct URCA process operates (which requires sufficiently large x, a condition 
necessarily satisfied at early times of the proto- neutron star), we have 

/ t0 = 1.13xl0 7 p?/ 3 T2 eV s- 1 , (14) 
C = 6.7x 10 21 pif g/(cms), (URCA) (15) 

where the £ expression is evaluated in the kq^> uo limit. We see from Fig. 2 that for 
T > MeV, the viscosity is less than 10 25 g/ (cm s) (depending on whether URCA or 
modified URCA processes operates), and thus the bulk viscosity cannot suppress 
the CFS instability of the bar-mode (see Lai 2001 for more details). 



IV EFFECT OF STAR-DISK COUPLING 

A newly formed rotating neutron star is often surrounded by a disk. As the bar 
grows in the NS, it will excite density waves in the disk, and therefore transfers an- 
gular momentum to the disk or from the disk. This will either enhance or suppress 
the bar-mode instability — this is simply another form of CFS instability. 

The angular momentum transfer is mainly through the so-called Lindblad res- 
onances, where 2(Q P — f^) = ±Qk (Q P = oJi/2 is the pattern rotation of the bar 
and Vtk is the rotation of the disk, assumed to be Keplerian). The torques can be 
calculated using the formalism of Goldreich & Tremaine (1979). At the inner Lind- 
blad resonance (ILR), Q k = 2Q P , the driving rate of the bar- mode due to star-disk 
coupling is (Lai 2001) 



6tt 2 n 2 p S(r ILR ) 
5 (to, - 2Q P ) (m/riy 



7ILR —J^ oo \ ~F~aH~I n?\ ? ( 16 ) 



where the negative sign in the front implies that the torque tends to damp the 
mode, Q s is the rotation rate of the star, R e is the equatorial radius, and E(riLR) 
is the surface density of the disk evaluated at ILR. Similarly, at the outer Lindblad 
resonance (OLR), we have 

98tt 2 g S(r OLR ) 
7oLR " 45 (fi s - 2Q P ) (M/Rl) ' 1 } 

where the positive front sign means that the torque drives the CFS instability. 
Clearly, the net effect of star-disk coupling on the mode depends on the relative im- 
portance of ILR and OLR. Compared with driving rate of the mode due to gravita- 
tional radiation, we find that star-disk coupling is important when £ > 10~ 6 M/i?g. 
So even a small amount of material outside the proto-NS may potentially affect the 
CFS instability of the bar-mode. 



V NONLINEAR EVOLUTION OF BAR-MODE AND 
GRAVITATIONAL WAVEFORMS 

In general, to determine the nonlinear evolution of the bar-mode requires 3d hy- 
drodynamical simulations including gravitational radiation reaction, and one needs 
to follow the system for a time much longer than the dynamical time of the star. 
This possesses a significant technical challenge (see recent attempt by Lindblom et 
al. 2000 on the r-mode evolution where an approximate ansatz for the radiation 
reaction is adopted). 
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FIGURE 3. The amplitudes of the GWs emitted by a secularly unstable NS, evolving from a 
Maclaurin spheroid toward a Dedekind ellipsoid. The two curves correspond to initial (3 = 0.24 
and 0.2 respectively (n = 1 is the polytropic index). The GW frequency sweeps downward. 



A Dedekind-like Bar 



One of the advantages of the m = 2 f-mode (as opposed the higher m f-modes 
and or r-modes) is that under certain idealized condition, namely for impressible 
fluid, there exists an exact solution for the nonlinear development of the bar-mode 
instability driven by gravitational radiation (Miller 1974; see Lai & Shapiro 1995 
and references therein; the latter also includes an approximate, compressible gen- 
eralization). 

The evolutionary sequence is as follows: We start with an axisymmetric Maclau- 
rin spheroid with f3 > 0.14 (secularly unstable). The bar grows and has a pattern 
angular frequency fl p (which is related to the mode frequency in the inertial frame 
by flp = Wi/2). Relative to the bar, there is also an internal rotation Q in which is 
larger than Q p . The important point to note is that although the mean rotation of 
the star Q s ~ Q p + Q in is near breakup, Q p can be much smaller (in fact, at the bi- 
furcation point (3 = 0.14, we have Q p = 0). As the amplitude of the bar continues to 
grow, the bar also gradually slows down. Eventually we reach a configuration with 
zero pattern speed. This is the Dedekind ellipsoid, basically stationary "football" 
with a fixed figure in space but with a lot internal rotation fij n . 

The gravitational wave (GW) emitted during such a quasi-equilibrium secular 
evolution is quite interesting. Figure 3 shows the GW amplitude as a function of 
GW frequency. The GW is quasi-periodic. Initially, we have an axisymmetric star, 
so h — 0. Then the GW amplitude increases as the bar grows. In the meantime, 
the bar slows down and the GW frequency (= 2Q P ) decreases. So eventually h 
decreases. Thus we have a non-monotonic GW amplitude evolution, with the GW 
frequency sweeps downward from a few hundred Hertz toward zero. The timescale 



FIGURE 4. The amplitudes of the GWs emitted by Jacobi-like bars. The two curves correspond 
to final (3 = 0.24 and 0.12 respectively. The GW frequency sweeps upward. 
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of the evolution is of order seconds to minutes, and the characteristic number of 
cycles of GWs is of order 10 4 . 

So far we have assumed the star begins its evolution from an axisymmetric 
state. Recent numerical simulations (e.g., New et al. 2000,; Brown 2000; Shi- 
bata et al. 2000), however, indicate that at the end of the dynamical evolution, the 
star may be already elongated rather than axisymmetric. So one can ask about 
the long-term evolution of the bar and the emitted GWs. Here the ellipsoid model 
can also provide qualitative answers. There are two possibilities (see Lai & Shapiro 
1995 for details): (1) If the bar has Q in > Q p , namely if the bar is Dedekind-like, 
then the evolution and waveform discussed above should also apply except that we 
need to cut off the initial growth phase of the bar. (2) Another possibility occurs 
when Q p > Q in , which we discuss in the next subsection. 

B Jacobi-like Bar 

If the bar has internal rotation (relative to the bar figure) less than the pattern 
rotation, i.e., Qi n < Q p , the bar is called Jacobi-like, and the evolution is quite 
different. In this case, gravitational radiation reaction tends to decrease the am- 
plitude of the bar, making the star less elongated. In the meantime, Q p increases 
because, even though J decreases, the moment of inertia decreases faster. So the 
GW frequency increases and eventually the star becomes axisymmetric. Here again 
we have a quasi-periodic GW signal (see Fig. 4), except that the frequency sweeps 
upwards. The timescale foe the evolution is of order a second, and the number of 
cycles is of order a few hundred. 

Note that at the end of the Jacobi-like evolution, the axisymmetric star may still 
be secularly unstable. If so the star will continue to evolve in the way described in 
Sec. V.A. * 



C Characteristic GW Amplitude 

To recapitulate, there are two types of rotating bars, and their evolution and 
emitted GWs are qualitatively different: (1) A Dedekind-like bar has large in- 
ternal rotation, and the resulting waveform sweeps downward in frequency; (2) A 
Jacobi-like bar has relatively small internal rotation, and the GW frequency sweeps 
upwards. Figure 5 shows the characteristic GW amplitudes for the two types of bars 
compared with the sensitivity curves (h rms ) of LIGO I, II and III. The characteristic 
GW amplitude is given by (Lai & Shapiro 1995) 



h c = h 
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where D is the distance, | dN/d In f\ = \f 2 /f\ is the number of cycles of GW spent 
near frequency f, E = E(M 2 /R ) is the energy of the star, and Q p = Q p (M / R^) 1 ' 2 



is its angular pattern speed. For a broad band detector such as LIGO, the best 
signal-to- noise ratio will be obtained by matched filtering of the data, with S/N ~ 
h c /h rms . Note that unlike coalescing compact binaries, the phase evolution of the 
GW from the evolving bars can not be determined with the accuracy needed for 
matched filtering, so such a high S/N may not be achieved in practice. A new fast 
chirp transform technique (Jenet & Prince 2000) is promising for detecting such 
signals. 

The event rate of Type II SNe at distance of 30 Mpc is about 100 per year. So 
even if a small fraction of the NSs are formed rotating rapidly, the GW signals 
discussed here are promising for 1IGO. Of course, it should be emphasized that 
the waveforms discussed in this section are based on the exact solution for the 
idealized situation (incompressible fluid). Whether this idealized solution has any 
resemblance to reality remains to be seen by future studies. 



VI PARAMETRIZED WAVEFORMS 

Finally we give fitting formulae for the gravitational waveforms generated by a 
rapidly rotating NS as it evolves from an initial axisymmetric configuration to- 
ward a triaxial ellipsoid (Maclaurin spheroid =^ Dedekind ellipsoid) as discussed in 
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FIGURE 5. Comparison between the characteristic GW amplitudes h c emitted during the 
secular evolution of a nonaxisymmctric neutron star and the rms noise h Tms of LIGO I, II and III. 
The solid and dashed lines correspond to the two different cases in Fig. 3 and Fig. 4. 



Sec. V.A. We use units such that G — c — 1. 

The waveform (including the polarization) is given by eq. (3.6) of Lai & Shapiro 
(1995) (hereafter LS). Since the waveform is quasi-periodic, we will give fitting 
formulae for the wave amplitude h (Eq. [3.7] of LS) and the quantity (dN/dln f) 
(eq. [3.8] of LS; related to the frequency sweeping rate), from which the waveform 
h+(t) and h x (t) can be easily generated in a straightforward manner. 

1. Wave Amplitude: The waveform is parametrized by three numbers: f max 
is the maximum wave frequency in Hertz, M 14 = M/(1.4M ) is the NS mass in 
units of 1.4M , R w = i?/(10km) is the NS radius in units of 10 km. (Of course, 
the distance D enters the expression trivially.) It is convenient to express the 
dependence of h on t through / (the wave frequency), with f(t) to be determined 
later. A good fitting formula is 
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where 



A = f (W1756) 2 - 7 , for f max < 400 Hz; 
" ~ 1 (W1525) 3 , for f max > 400 Hz; 

Note that if we want real numbers, we have 

M 2 99 9 , /30Mpc 

— = 4.619 x lO^M^i?- 1 [-^- 

2. Number of GW cycles: The fitting formula is 



with f max = f m axR% 2 M 1 l /2 . (20) 
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Notes to the fitting formulae: 

Note (i): Using the above equations, we obtain the characteristic amplitude: 
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which agrees with Eq. (3.12) of LS to within 10% [Note that in Eq. (3.12) of LS, the 

factor f 1 / 2 should be replaced by (R^qMi^ 2 f) 1 ^ 2 , similar to the above expression.] 
Note (ii): The accuracy of these fitting formulae (as compared to the numerical 
results shown in LS) is typically within 10%. When / is very close to f ma x, the 
error in the fitting can be as large as 30%. 



Note (iii): The frequency evolution fit) is obtained by integrating the equation 
f 2 /f = — \dN/d\nf\ (note that the frequency sweeps from f max to zero). For 
example, we can choose t — at / = 0-9f max . (Note that one should not choose 
t = at / = fmax as the time would diverge — the actual evolution near f ma x 
depends on the initial perturbations). Once fit) is obtained, the waveform can be 
calculated as (cf. Eq. [3.6] of LS): 

1 4- cos 2 9 

K = h[f(t); f max , M, R] cos $(t)^ , (24) 

hx = h[f(t); f max , M, R] sin cos 9, (25) 

where 9 is the angle between the rotation axis of the star and the line of sight from 
the earth, and $(t) = 2n f f(t)dt is the phase of the gravitational wave. 

Note (iv): f max typically ranges from 100 Hz to 1000 Hz (see Fig. 5 of LS); M 1A 
and R 10 are of order unity for realistic neutron stars. 
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